Aging dynamics of quantum spin glasses of rotors 
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We study the long time dynamics of quantum spin glasses of rotors using the non-equilibrium 
Schwinger-Keldysh formalism. These models are known to have a quantum phase transition from 
a paramagnetic to a spin glass phase, which we approach by looking at the divergence of the 
spin relaxation rate at the transition point. In the aging regime, we determine the dynamical 
equations governing the time evolution of the spin response and correlation functions, and show 
that all terms in the equations that arise solely from quantum effects are irrelevant at long times 
under time reparametrization group (RpG) transformations. At long times, quantum effects enter 
only through the renormalization of the parameters in the dynamical equations for the classical 
counterpart of the rotor model. Consequently, quantum effects only modify the out of equilibrium 
fluctuation dissipation relation (OEFDR) , i.e. the ratio X between the temperature and the effective 
temperature, but not the form of the classical OEFDR. 



PACS: 75.10.Nr, 05.30.-d, 75.10.Jm, 75.lQ.Hk 
I. INTRODUCTION 

In recent years the focus of studies of glassy systems 
has moved from equilibrium to non-equilibrium proper- 
ties, reflecting the understanding that glassiness is an 
intrinsically non-equilibrium phenomenon. The effect of 
quantum mechanics on spin glasses has also been a topic 
of interest. On the experimental side, there have been 
several studies of spin systems which are in the vicin- 
ity of a cwajjitum phase transition from a spin-glass to a 
spin-fluido □ including tie dipolar, transverse-field Ising 
magnet Li Ho^: Yi-x On the theoretical side, studies 
have mainly focussed oijit|t|h<; Ising model in a transverse 
field or quantum rotorsudlj 

There has also been much recent vrorkjOpi. h ft q uantum- 
classical connection in spin glasses jl3'tjt3'E3Eil such as 
the derivation of quantum TAP equations ,1111 and tha 
large N solution of SU{N) models of spin glasses.Ell 
Quantum fiuctuations have been shown to be responsible 
for redefining the boundary between the spin glass and 
the quantum disordered phases in systems such as the ro- 
tors or the transverse field Ising model. In systems such 
as the quantum extension of the spherical p-spin-glass, 
quantum fluctuations arc responsible for a-erossover from 
a second to a first order phase transition.111 

Another system in which quantum effects and glassi- 
ness coexist is in the stripe glasses found in doped Mott 
insulators (e.g. Ref. p^ . 

In principle, quantum effects could also alter the non- 
equilibrium dynamics in the glassy regime at very low 
temperatures. Whilst the properties near the T = 
quantum critical point of these S5p±ems can be under- 
stood within a static formulationjlj studying the non- 
equilibrium effects in the glassy system requires a quan- 
tum dynamical approach. 

Aging effects in classical spin glasses and other glassy 
systems have been investigated considerably in recent 



years. The out of equilibrium nature of the system dis- 
plays itself through a persistent dependence of suscep- 
tibilities on the waiting time,f-t|he time since the sys- 
tem entered the glassy phase. E3 The correlation func- 
tions in sucb-systems often display quite general scal- 
ing features,L3 and mayj-also have an invariance under 
time reparametrizations.cJ It is a natural question to ask 
whether quantum spin glasses have a similar aging be- 
haviour to their classical counterparts. 

Understanding aging dynamics in quantum sys- 
tems requires formulating the problem within a non- 
equilibrium closed time path (OTP) formalism, such as 
the Schwinger-Keldysh approach. This dynamical ap- 
proach is particularly attractive for studying disordered 
systems, since it eliminates the need to use replicas in 
carrying out the average over quenched disorder. This 
is because the generating functional in the CTP formal- 
ism is independent of the disorder realization. Several 
different disordered quantum systems have been studied 
recently using this approacL, e.g. disordered and inter- 
acting electronic systemsJHj'Ej the iafinite-range quantum 
p-spin glass in the spherical limit ,L3 and a semionic rep- 
resentation for quantum spin systems.EZ 

In this paper, using the Schwinger-Keldysh dynamical 
approach, we investigate the aging behaviour of a spin 
glass of M component quantum rotors that have infinte 
range interactions and on site self interactions with a 
coefficient u. The model of quantum rotors is consider- 
ably simpler than that of true quantum Heisenberg spins 
present in any isotropic antiferromagnet or systems like 
the doped cuprates; the different components of the rotor 
variables all commute with each other, unlike the quan- 
tum spins. As a consequence, the path- integral written 
in the rotor variables has an action which contains no 
Berry phases and is purely real. 

The quantum rotors are particularly interesting in elu- 
cidating the role of quantum mechanics in the long time 
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dynamics of the glassy phase for a number of reasons. 
They provide an example of a quantum system with con- 
tinuous replica symmetry breaking (RSB), in contrast to 
the one-step RSB quantum p-spin model with a spheri- 
cal constraini- which was recently studied by Cugiiandolo 
and Lozano.Ej Another feature of the rotor system is that 
it contains an expansion parameter, u, which we can use 
to organize a perturbative expansion. When we consider 
the model to 0(u), we get equations that look very sim- 
ilar to those obtained in the case of the p = 2 spherical 
spin model, which has no RSB. However, when we ex- 
pand to 0{u^) we end up with terms that are similar to 
those from the p — 4 spherical model and also some terms 
that do not arise in p-spin models. Therefore, we can in- 
vestigate very generally how each term contributes to the 
long time dynamics of the model, and how the quantum 
effects enter in the problem through each of these terms. 

We show that the terms in the dynamical equations 
that appear as a consequence of quantum mechanics are 
irrelevant at long times. The precise sense in which the 
terms are irrelevant was defined in Ref. using the 
language of time reparametrization invariance and the 
reparametrization group (RpG) of time transformations. 
Hence, at very long times, the dynamics of the quantum 
rotors is completely determined by a renormalized clas- 
sical version of the model. In this time regime, quantum 
mechanics enters the problem only through a renormal- 
ization of the coefficients in the dynamical equations, and 
there is a complete correspondence between the classical 
and quantum versions of the model in the aging regime. 
In particular, this implies that the response and corre- 
lation functions of the quantum rotors are related by a 
modified fluctuation-dissipation relation, with an effec- 
tive temperature Toff > T, just as found in the classical 
spin glass models. Our results within the UpG frame- 
work extend the connection between the aging regimes 
of quantum and classical systems, as found in the quan- 
tum p-spin modelMI with one-step RSB, to a wider and 
more general class of quantum glassy models. 

The paper is organized as follows. In section |l|, we de- 
scribe the model that we study and the CTP formalism 
that we use. In section [II we obtain the saddle point dy- 



namical equations for the model and develop the pertur- 
bation theory in u that we need to obtain these equations. 
In section IV we study the solutions to these equations in 
the paramagnetic phase. In section ^ we investigate the 
solutions to the dynamical equations in the spin glass 
phase, and show that a selected number of terms can 
be made invariant under reparametrizations of the time 
coordinate in the aging regime. We analyze the terms 
stemming from the quantum dynamics imposed on the 
system that are not present in the classical model, and 
show that they become irrelevant at long times under 
the time reparametrization transformations. Finally, in 
section ^ we discuss differences between classical and 
quantum spin glasses. 



II. THE MODEL 

We consider the model of a glass of quantum rotors in- 
troduced in Refs. with long range interactions on a 
d-dimensional lattice with iV sites. (We let d and N go to 
infinity later) . An important point to note is that the M 
components of angular momentum of the rotors commute 
on the same site, unlike Heisenberg spins, for which there 
are non-trivial commutation relations between the com- 
ponents. This simplification allows us to write a path 
integral in which there are no Berry phase terms, and 
makes the rotors much easier to treat analytically. 

To derive the dynamical equations ■ , \vjs , use a CTP 
(Schwinger-Keldysh) formalism. yoy'BEl The CTP 
approach provides a way to study the non-equilibrium 
response of a system. The price that has to be paid for 
this is the introduction of a second component to the 
system, with time flowing in the opposite direction. The 
non-equilibrium formulation has the advantage that it 
leads to a generating functional that is automatically nor- 
malized to unity. The property of normalization allows 
averages over disorder realizations in a way that bypasses 
using replicas. 

There are many possibilities for the ApM^*^ integra- 
tion contour in the complex time plane.EJ'Eil The contour 
C that we work with is illustrated in Fig. [l], which starts 
at t = 0, runs to t = cx), and then returns in the negative 
t direction to t = 0. The contours more usually used, 
for example in problems involving the calculation of non- 
linear. response, run from t = —oo to t = cxd and back 
againJ23i£j however after an infinite time, the system will 
have equilibrated, so it is no longer possible to study the 
non-equilibrium dynamics that we are interested in. 



FIG. 1. The contour C used in the formalism. 

For a system at equilibrium, the usual way to intro- 
duce finite temperature is through the density matrix 
with the form of the Gibbs-Boltzmann distribution. In 
a non-equilibrium situation, the system cannot be de- 
scribed by the Gibbs-Boltzmann distribution, so an al- 
ternative approach is required. The solution is to couple 
the system to a heat bathE3 (chosen here to be a set of 
independent harmonic oscillators) and then allow the sys- 
tem to reach a constant temperature. A detailed account 
of integrating the bath variables to obtain an effective ac- 
tion in the p-spin model is given in Ref. Ua, and we follow 
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their approach in our study of the quantum rotors. 



where the noise and dissipative kernels, v and rj respec- 
tively, are given by 



A. Effective action for the rotors and bath system 

The action for the system of rotors interacting with 
the bath takes the form 



S — 5frco + S'int + Sdis + St 



(1) 



where S'fice is the free action, ^int describes the self- 
interaction of the rotors, Sdis contains the spin exchange 
interactions which introduce disorder and frustration to 
the model, and St describes the interaction with an ex- 
ternal heat bath. The Lagrangian for the free rotors is 



—m^ Sf^ 173'' S^^ 



(2) 



The indices on the spin variable Sf^ refer to the Keldysh 
contour branch (a), site (i), and spin component (/x); 
sums over repeated Keldysh and spin component indices 
are always implied unless explicitly noted otherwise. The 
CTi 2,3 are the standard Pauli sigma matrices. Note that 
the S^ terms enter with an opposite sign to the S^ terms 
because the direction of time integration is reversed on 
the second part of the contour. The first term in the 
Lagrangian Chcc is a kinetic energy term {1/g is the mo- 
ment of inertia), and the second is a potential energy 
term (m^ also acts as the parameter that tunes between 
the paramagnetic and spin glass phases) - the free La- 
grangian may be thought of as M-component harmonic 
oscillators on a lattice, each at site i. 
The self-interactions are described by 



-Cint — 



(3) 



The strength of the self-interaction is u and this term is 
necessary to ensure the stability of the system, since oth- 
erwise the Lagrangian, would describe randomly coupled 
harmonic oscillators .Ej 

The term in the Lagrangian that contains quenched 
disorder is 



na . 



S 



(4) 



Thermal effects are introduced by placing the system in 
contact with a heat batii, using the method introduced by 
Feynman and Vernon.E3 The contribution of this coupling 
to the action after integrating out heat bath variables is 

St= - I dhdh J2 Vih - t2) Slit,) [at + laf] Slih) 

i 

I J dhdt2 - h) Sf^ih) - af] S^ih), 

i 

(5) 
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y{t)= J duj I (uj) coth I— ^ ] cos{ujt), (6) 



(7) 



T]{t)= -e{t) I duj I{uj) s\n{ujt), 
Jo 

and I{uj) is the spectral density of the bath 



lie) 



E 



SiuJ - LOn): 



(8) 



where Nb is the number of oscillators in the bath, w„ 
is the natural frequency of the n*^ oscillator, Af„ is its 
mass, and Cn is its coupling to the system. Here we only 
consider the case of ohmic dissipation 



/(cj) = —uje 



for < A, 



where 7 plays the role of a frictioiiL.cQefRcient (to relate 
to the notation in previous worksjl3'E3 7 = M70). 

Properties to note are that both kernels are purely real, 
I'it) = i^i—t), and both kernels decay rapidly at large 
time differences. Note that the quantum fluctuation- 
dissipation theorem (QFDT) holds for the bath variables 
when there is ohmic dissipation. The QFDT relating the 
heat bath variables ialj 



, . 1 fdco' 
r]iuj) = - hm 



2tt lu ■ 



■ tanh I — - — 1 nvyLo ). 



(9) 



B. Disorder average and definition of correlation 
functions 

Having described all of the terms in the action, we 
now move on to performing the average over disorder 
realizations and defining the correlation and response in 
terms of the CTP two-point correlators. 

The closed time path generating functional is 

Z = 1 \DSt^\ e^^ . (10) 
We perform an average over the quenched disorder 



Z= \ \DJ\ViJ,,) [DSUe 



where 



ViJ^,) = 



/ N 



V27rJ2 



exp 



NJl 



2J2 



(11) 



(12) 
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for a long range interaction. The form Eq. (|T^ for the 
disorder probabihty distribution comes from the assump- 
tion that the disorder is Gaussian distributed and that 
Jfj — J'^/N . We can perform the integration over disor- 
der without rephcas because of the normahzation prop- 
erty of the generating functional in the CTP formulation, 
since the disorder and the initial conditions are uncorre- 
lated: 



I [DJ] V{.h,) e^/'^'S.^. < = 

(13) 

The dynamical equations are written in terms of the cor- 
relation and response, which are defined below. We use 
an overbar to indicate an average over realizations of 
disorder and angular brackets (. . .) to indicate an aver- 
age with respect to the action. The correlation, C(ti, t2), 
and response, R{ti,t2), are 

1 



a.Mtuh) - ^{suh)sut2) + sfJh)sUh)), (14) 



and Rij,fiv{ti,t2) = 5{S}f^{ti)) / 5hj^{t2), which in linear 
response theory may be written as 



R 



.,,Mti,h) = ^{Sl{h)[S}Jh) - S],ih)]) . (15) 



An alternative approach to the dynamics is to intrp-. 
duce a matrix propagator as in the Keldysh formalism,c2l 
which takes the form for bosonic fields (e.g. Ref. |3^ ) 

Gll^,{h,t2) = -i{T [Slih)Sj,{h)]) , (16) 

Gll^,{h,t2) ^ -i{Sl{h)Slih)) , (17) 

G^l^AhM) = ~i{Slih)Slit2)) , (18) 

Gll,Ati,t2)-~^(f[Sl{h)Slit2)]), (19) 



where T and T are the operators for time ordering and 
anti-time ordering respectively. Under the transforma- 
tion in Keldysh space G ^ G = LGL\ where 



a/2 



(1 - 10-2), 



the propagator takes the form 
G = LGL'' = 



G^ 



(20) 



(21) 



where 



and the results follow from the definitions of the propa- 
gators in terms of time ordered products. The definitions 
above imply the following properties for the propagators 



Gft,i,^{t2,tl) — Gfji^^^{ti,t2), 



Gji,zy^(^2, ^l) 



uih,t2). 



(25) 
(26) 



The retarded and advanced propagators are purely real 
and the Keldysh propagator is purely imaginary. Under 
the transformation, L, the spins are also rotated, 5* 
S = LS, so for 



S 



we have 



V2 \S' + S' 



(27) 



(28) 



The S and S variables in the quantum CTP formulation 
naturally become the usual variables within the Martin- 
Siggia-Rose (MSR) formalism^ when the classical limit 
is taken. The relation between the Keldysh and retarded 
propagators and the correlation and response is 



G'Hh,t2) = -2iG{hM): 
G'^ih,t2) = -hRih,t2). 



(29) 
(30) 



We will work with the correlation and response rather 
than the propagators in the dynamical equations. How- 
ever, it is easier to use the matrix propagator in the Feyn- 
man rules for including interactions, rather than using a 
diagrammatic technique where response and correlation 
are treated differently.E3 

In anticipation of the saddle-point evaluation for the 
disorder term, we perform a Hubbard-Stratonovich trans- 
formation to decouple the four spin term generated by the 
disorder average. Introducing the Hubbard-Stratonovich 
field Qij^f^„{ti,t2), the effective action can be rewritten 
as 



Z 



Z[Q] = J [D^]ei^^f-=+i^^Q+i-5.nt-HiST^ 



where 



(32) 



^ftcc ^ JdtJ2 ^^f^ (*) Si (t) , (33) 

i 

Sq = j dt,dt2Y,Qll,{t,M) ~Sl{ti) ~SUt2) . (34) 

i 

Next we introduce a modified version of the Hubbard- 
Stratonovich field, Qf^^ = Qf^^^ erf, thus 



G^ 


= G" 


-Gi2 


= G2i 


-G2^ 


(22) 


Z 


G^ 


= G" 


-G^i 


= Gi2 


-G22, 


(23) 




G^ 


= G" 


+ G22 


= G2i 


+ G12, 


(24) 





1) Q°„.(tl.*2) 



(35) 
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where 



The f term has the structure 

f afc _ r ( 



(36) 



where 



(37) 



(38) 



Sj(ii,i2) 



2J2 



(Tl G{ti,t2) CTl 



G« 




(42) 



(tl,t2 



The notation of a matrix with (ti, f2) following it is used 
to indicate that all of the propagators in the matrix have 
that as their argument. At the mean field level the so- 
lution is homogeneous in the site index, so we drop it. 
In the absence of a magnetic field, we also drop the spin 
indicies, i.e. Gp,,y — G 6^,^- 



B. Diagrammatic perturbation in u 



III. DYNAMICAL EQUATIONS 

Using propagators, the self-consistent mean field equa- 
tions can be represented in the form 



Gn. 



(39) 



where the S terms are the self energies from the disor- 
der, the interaction with the thermal bath and from self- 
interaction respectively. The general strategy that we 
will adopt in obtaining the dynamical equations for this 
model is to perform a saddle point evaluation of Ej and 
then a perturbation expansion in u. We obtain a solution 
to 0{u) in the paramagnetic phase, in analogy with t 
replica symmetric solution in the equilibrium problem, 
whilst in the spin glass phase we need to consider terms 
to 0{v?) in the interaction, which are the terms found to 
contribute to RSB. 



A. Saddle point evaluation of Ej 

At the saddle point the functional derivative with re- 
spect to the Hubbard-Stratonovich field is zero 

5Z 

The variation leads to the following equations 


+E('5S^(^i)^i(*2)). (40) 

i 

Rearranging gives 



/ hb' nb'a' I, , \ a' a 



Qt^^ihM)) - -^- of G^';'(^2,<i) (41) 



By substituting the saddle solution into Eq. (pq), one 
obtains the self-energy contribution coming from the dis- 
order: 



The interaction terms in the action may be treated in 
a diagrammatic way by noting that in the S basis, the 
interaction term is written as 



dt St^it)SUt)S^t)a1''Stit). (43) 



The Feynman rules (see Fig. g) for calculating E„ term 
by term in u are that the propagator (solid line) is 
iG'J;^'l{t,t') and the interaction propagator (dashed line 
with a dot at one end) is i^af'5'^'''8p,u- The ai is inserted 
at the end with the dot. 



a |i 



b V 



iG;;^(t,t ) 
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FIG. 2. The Feynman rules for the pertubation theory in u. 

This leads to four contributions to the self energy 
to first order in u, as shown in Fig. |[ Note that the 
bottom two diagrams in Fig. |^ are proportional to M, 
the number of spin components of the rotors, since there 
is a trace around a closed loop, and that the trace in- 
cludes spin indicies. Note also that tr{Gai) — 0, since 
G^{t, t) = 0. The sum of these contributions is 

= f (1 + M) a,G''{t,t)S{t~t') , (44) 

where the superscript (1) indicates that the self energy 
is to first order in u. 

To second order in u there are six different one parti- 
cle irreducible diagram topologies to consider which con- 
tribute to the self energy and four different diagrams for 
each topology (the different topologies are shown in Fig. 

We treat each topology separately, starting with the 
propagator with two self interactions, which contributes 



"[3(G^)2 + 3(G«)2 + (G^)2]G^ 
[3(G^)2 + (G^)2]G« 



[3(G^)2 + (G-4)2]G-4~ 




(45) 



(t,t') 
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FIG. 3. Diagrams contributing to the self energy to first 
order in u. 

The next diagrams are the vertex correction to the 

closed loop that contributes to si^'' and the propagator 
with one self interaction, which contribute 



E(2f)(i,t') 



4h 



CTl 6{t - 1') 
diG^il i)G"it, i)G^{t, i) 



(46) 



The remainder of the self energies are expressible in 
terms of the expressions found in equations ( [45|) and 

(H): = MS(2'^), ^ Afl](2b)^ S(2e) ^ j^(26) 

and S^^^) = S*^^''). The diagram that plays the most im- 
portant role in the physics is (2a) because it depends on 
two times rather than only one. 



e 



2a 



2c 



2b 



2d 



2e 



2f 



By substituting the self-energies Sj, St and S„ that we 
determined above to order m^, and shifting notation to 
response and correlation, we obtain the dynamical equa- 
tions for the quantum rotor model 



1 d 



1 



S{ti - i2) = - -m' - u{l + M)C{hM) 

-2{1 + M)u^C{ti,ti) 



10 

-2J2 



X / dtR{ti,t)C{ti,t)jR{ti,t2) 

dtT]{ti ~ t)R{t,t2) 

dtR{ti,t)R{t,t2) 

C{ti,tfR{ti,t)R{t,t2 



o roc 
(l + Af)M^ / dt 

2 Jo 



dtR{ti,tfR{t,t2), (48) 



and 



1 92 



-m^ -u{l + M)C{ti,ti) 



2g dtj 2 ■ 

2{1 + M)u^C{ti,ti) I dtR{ti,t)C{ti,t)}C{ti,t2) 



/o 

-2J- 
1 



/ dt7j{h-t)C{t,t2)- - / 

Jo ^ Jo 



dtviti - t)R{t2,t) 





/"OO 

2 ' [C(tl,t)i?(t2,<) +i?(tl,<)C(t,t2)] 



dtC{ti,tfR{t2,t) 



dtC{ti,tfR{ti,t)C{tM) 



^ -10 



dt [R{ti,t)^C{t,t2 
SR{ti,tfR{t2,t)C{ti,t)] . 



(49) 



The terms that arise in this expansion look partly like a 
p = A spherical spin glass (the C^R term in the corre- 
lation equation), but there are also terms (multiplied by 
h^) which do not appear in the classical spherical model. 



FIG. 4. Diagrams contributing to the self energy to second 
order in u. 



C. The dynamical equations 

3) we obtain the dy- 



From the self consistent Eq 
namical equation 



1 = (Go^i - Sj - St - S„) G 



(47) 



IV. PARAMAGNETIC PHASE 

In the paramagnetic phase the correlation and response 
are time translation invariant (TTI) after the initial tran- 
sients have died out. We find a solution to the dynamical 
equations ( ^8| ) and ( ^ ) to 0{u) in this phase. With the 
assumption of TTI correlators, let t = ti — i2, and then 
let ti 00. Define Goc = hmt^^oo C{ti,ti), which is the 
equilibrium limit of the equal time correlation (we also 
use Goo in the spin glass phase; as a one time quantity 
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the equal time correlation has a limit as ti oo). TTI 
allows us to solve the problem in the paramagnetic phase 
by performing a Fourier transformation: 



1 



1 



l— + -m^+u{l + M)C^ 

duo' R{uj')r]{u;') 
27r uj — uj' + iO 
duj' R{lj')R{lj') 
, 27r a; — cij' + iO ' 



(50) 



and 



dJ r]{ijj')C{uo') 
J 27r — + iO 
°° dm' v(uj')R{-uj') 
„ 27r uj — uj' + iQ 



-—I 
2 



duj' 
"2^ 



C{io')R{-uj') + R{uj')C{lo') 
w — + iO 



(51) 



The following properties hold for the correlation and 
reponse and their Fourier transforms: C{t) is invariant 
under time reversal and is real, hence C((jj) = C*(— a;), 
C*(w) — C{iu) and C{uj) is real. The response is real, 
hence i?(w) + i?(-w) = 2 Re[R{uj)] and i?(w) - i?(-w) = 
2iIm[i?(Li;)]. In addition, the noise kernel i^(w) is real 
and has the property iy{Lo) ~ i>{—lo). The real part of 
the equation (p^) combined with the QFDT for the heat 
bath variablesTfl) gives 



Im[i?(c^)] = itanh('^ 



(52) 



Hence the correlation and response are related by the 
QFDT if the same is true for the heat bath variables. In 
this regime, we can solve explicity for the response, and 
hence the correlation through the QFDT (|5^). We use 
these solutions to investigate the self-consistency condi- 
tion in the paramagnetic regime and investigate critical 
slowing down as discussed by Sompolinsky.E3 



A. Response, correlation, and phase boundary 



(note that ?7(w) has both real and imaginary parts). In 
the limit that w — > 0, the real part of the kernel rj is 
proportional to the frequency cutoff. A, so to remove this, 
define — — 2rj{0). Then at zero frequency at the 
critical point, 



+ u{l + M)Coo = 2^2 7, 

and the imaginary part of the response is 
Im[?7(w)] 



(54) 



lni[R{uj)] 



4J2 



1 



-^^2 Vlni[^7(t^)]2 +8J2 . (55) 



To get the self-consistency condition for Coo at the criti- 
cal point, we use equation (p5[) and the QFDT 







/ — coth 


J-oo 27r \ 





Using 



lm[r]{iu)] ^ _^e-l'^l/^6l(A- 



(56) 



(57) 



and noting that the second integral vanishes, since the 
integrand is odd, we get the result 



0(A2 



(58) 



The quantity Coo fixes the equal-time correlation, and 
provides a constraint for the rotor size in this .model, 
since we do not impose a spherical constraint .ESeZI In the 
spin glass phase C{t,t) should relax to an equilibrium 
value, q, the Edwards- Anderson order parameterH This 
is because it is a one time quantity ~ the same is not true 
for two time quantities. The quantity m\ acts to "tune" 
the system through the spin glass transition. Using equa- 
tion ( p^ and separating into a piece that depends on 
the cutoff and a temperature dependent piece, we have 
an expression analagous to the one found for rc{T) in 
the equilibrium study of quantum rotorsE^I 



u{l + M)7 



(59) 



Assuming that the integrands in the resonant integrals 
in equations ( pO| ) and ( |5l| ) fall off sufficiently quickly at 
large w, then they may be rewritten and solved for R{lo) 
to give 

i?(c.) = -^-^V^(c.)2-8J2, (53) 

where 

^2 2 

x{uj) = — + -m^ + it(l + M)Coo. - 7?(t^), 
2o 2 



B. Critical slowing down 

The dynamic transition from the paramagnetic to the 
spin glass-^phase was investigated by Sompolinsky and 
ZippeliusEj However, in their treatment they took the 
infinite time limit before, rather than after the N ^ oo 
limit, the opposite order of limits than that considered 
here. This leads to finite energy barriers between traps 
and hence an infinite hierachy of time scales and an er- 
godic solution. In the situtation we consider the system 
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is confined to a single ergodic comppnent and we are 
studying the relaxation within a trapO 

To connect with the work of Sompolinsky and Zip- 
pelius, define 



which leads to r(a;) ~ 1/w at small frequencies and 



(60) 



dri_ _ uj_ 

which is analagous to their result. 



C. Alternative approach — Schwinger-Keldysh 
Landau Theory 



In previous work on quantum rotor systems)130 the 
approach taken has been to obtain a replicated field the- 
ory in the Haibbard-Stratonovich Q fields. In this and re- 
lated work,E£l we have obtained directly dynamical equa- 
tions for spin correlation functions. In the following we 
will briefly show how to set up a Landau theory for the 
fields Q within a completely dynamical approach, and 
connect to both the replica results of Ref. |l2|,^ and the 
results obtained in the preceeding sections. 

In the dynamical Schwinger-Keldysh approach, if the 
system does reach equilibrium with the bath (as is the 
case in the paramagnetic phase), then the temperature 
of the heat bath can be introduced through the use of 
a density matrix rather than with the heat bath kernels 
referred to in the previous section. 

Just as in section |l|, integrating over disorder and in- 
troducing a Hermitian Hubbard-Stratonovich field leads 
to the following action 



where Zq is the single site generating functional 



(62) 



So= Sfrcc + Snl (63) 



where Sni includes the self-interactions and interaction 
with the bath. We can expand Zs[Q] by obtaining the 
perturbative vertices in powers of the Q fields, noting 
that the mean-field solution is isotropic in space. Using 
standard diagrammatic methods and in simplified nota- 
tion, the effective action up to the cubic order is 



Tr[(GogGo)ai(GogGo)] 
Tr(GoQGo) 



-lU 



dioi du)2 dio^ 
"27r '2tt "27r" 
dwi duJ2 duj3, 
'2n "27r" 
Tr(aiGoQGo) + . . . . 



-lU 



(64) 



where . . . means higher order terms. It is important to 
point out that terms with coefficient u are important to 
get stable saddle point solutions. We study the saddle 
point solutions in the absence of spatial fluctuations. 

Consider the following saddle point solution which is 
0(M) and TTI: 

Qf,iuji,uj2) = 2nS{uji+LJ2) S^, Q°'(wi) , (65) 



Qk Qb. 
Qa 



(66) 



This Ansatz has the same structure in—frequency and 
0{M) indices as in the replica approach. However, here 
we are dealing with a 2 x 2 matrix in Keldysh space in 
contrast to the ti x n ( n — > limit ) replica matrix. This 
leads to the saddle point action 

Ssp ^ f (QrGr + QaG^) 



duj 
2^ 



-^u{l + M) J ^ [QniGir + QAiG'^n 

X / ^ [{QrG% + QaG'^)G% + QkG%G'^] 
+ ... . (67) 

By using the analytic properties of retarded and ad- 
vanced Green functions, we can remove the quadratic 
term by a shift Q{oj) Q (oj) — CT (uj) . In the low energy 



r, Gl 



K 



limit, using free propagators: G^ ~ G^ 
0, we get: 

5sp = J^{^'-r){Qn + QA) + ^J^[Q% + Ql] 
- iu{\ + M)J^ [Qr + Qa]J^Qk + ... . (68) 

Solution in the Paramagnetic phase 

Minimizing Equation ( |6^ ) with respect to Qr and Qa, 
gives 

= 0, (69) 

/I 
-^Qk = 0- (70) 
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The two equations are complex conjugates of each other, 
therefore there is only one independent equation. Solving 
for Qr gives 



(71) 



r^r-u{l + M) I — coth i^), (72) 

(73) 



X"{uj) - sgn(a;)(tj2 _ f)i/2gi(^ _ 



The critical point is located at f = where r = rdT) — 
u{l + M) J ^ Lu coth These results are the same as 
those obtained by the replica approach in Ref. HSl 



V. GLASSY PHASE 

In the spin glass phase the correlation and response 
are no longer time translation invariant. Instead, there 
is an aging piece that remains a function of two times. 
Below we investigate the dynamical behavior of the cor- 
relation and response functions for the quantum rotors 
in the glassy phase. 



A. Weak ergodicity breaking and weak long term 
memory 

Weak ergodicity breaking (WEB) and weak long term 
memory (WLTM) are phenomena that have been ob- 
served iiL.|Lhe_aolutions to classical mean field spin glass 
models.LjCjO The numeidcal results for the quantum 
version of the p-spin modelij are in agreement with those 
predicted from these scenarios, and we also assume WEB 
and WLTM in our solution of the quantum rotor model. 

WEB can be summarized as follows: the correlation 
function behaves as the sum of a TTI piece and a piece 
which depends on two times 



C{ti,t2) = CsT[tl~t2) + CAG{tlM) 



(74) 



The stationary (TTI) piece decays to zero as ii — ^2 oo, 
and the aging piece satisfies CyiG(i, t) = q, and 



lim CAaiti^h) 



0, 



(75) 



so that both times are important in the decay of the aging 
piece of the correlation. We separate the response in a 
similar way to the correlation 



R{tiM) ^ RsT{tl - t2) + RAGitlM) , 



(76) 



The stationary piece decays to zero as ti— 12 oo, whilst 
the decay of the aging piece depends on both times. The 
assumption of WLTM applies to the integral of the re- 
sponse function: the integral of the response over any 
finite time interval vanishes, i.e. 

lim [ dt2R{ti,t2) = , (77) 



for fixed t, however the integral over an interval that 
grows with time is finite 

lim / dt2R{ti,t2) ^0 . (78) 



B. Modified QFDT 

Treating the spin glass regime in a similar way to the 
paramagnetic phase, and looking in the limit ti ^ oo, 
it is possible to get an equation that includes the TTI 
pieces, and also some static pieces due to a non-zero 
Edwards- Anderson order parameter. 

In the spin glass phase, we need to take into accpuHt 
the aging that has occurred up until the waiting timeJl3'tII 
To do this we split R and C into stationary and aging 
pieces before letting ii — > oo and then Fourier transform- 
ing. The relation we find between the stationary parts of 
correlation and response is 

Im[i?sT(w)] = ^^^^^ (^~^) ^■S'r(^) 

+ ^^q iXco + {RC)oo) S{lu) , (79) 
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where 



Xoo = lim / dtRsT{ti-t), 
Jo 



(80) 



and {RC)oo is (noting that ii ^ t2 ^ |ti — ^2! in the 
evaluation of the integral) 



(i?C7)oo=lim / dtRAG{ti,t)[CAG{t.ti) + CAG{ti,t)) 



(81) 



Equation ([79|) shows how the QFDT is modified in the 
spin glass phase for the rotor model. There is the usual 
piece that is present in the paramagnetic phase, and then 
a delta function at zero frequency that grows with de- 
creasing temperature. 



Modified QFDT within the dynamical Landau theory 

As in the treatment above, in the spin glass state, the 
Keldysh component Qk acquires a non-trivial 5 function 
part due to a nonzero Edwards- Anderson order parame- 
ter 

iQK - iQk^ + 2TTq5{oj). (82) 
Substituting the ansatz into Equation ((7(1) gives 
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iv^-r + Ql- iu{l + M)J^ g^^s -{u + Mv)q = 0, 

(83) 

and solving for Q r and q leads to 



1 



q = 



u{l + M) 



{rc{T) - r). 



(84) 
(85) 



Again, these results are the same as those obtained by 
the replica approach in Ref. ^ 



C. Aging regime 

Aging behaviour has been shown to be a feature in 
the long tirne. dynamics of many glassy systems, both in 
experimcntfaS and theoretical models (see Ref. 22 for 
a review). In the case of spin glasses, the major models 
that have been studied in the aging-jregime are the Sher- 
rington_KirkpatrickE3 (SK) modelJ13 the spherical p-spin 
model, ^ and the quantum version of the spherical p-spin 
model.Ej In our work we consider the quantum version 
of the rotor model, which in the absence of a magnetic 
field, is similar to the soft spin SK model with quan- 
tum dynamics rather than Langevin dynamics, having 
M components instead of only one. 

One of the features of equilibrium solutions of mean 
field spin glass models is replica symmetry breaking, in- 
troduced by ParisieZl in the context of the SK model. 
Cugliandolo and Kurchanc3 started from a set of equa- 
tions for the long time dynamics of the SK model and 
showed that the correlations have an ultramctric struc- 
ture analagous to that found in the Parisi RSB solution 
of the SK model. This ultrametric structure in the long 
time dynamics can be obtained in the model here only by 
including terms to O(w^) in the dynamical equations - 
there is a direct Malogy with the equilibrium solution of 
quantum rotors,t2l where RSB occurs when 0{v?) terms 
are included in the equilibrium solutions. (The approach 
of Read, Sachdev, and Ye was slightly different from ours, 
since they expressed all terms in the action in terms of 
the Hubbard-Stratonovich field Q°'^ , rather than working 
directly with the spin correlations, as we do here). 

An important feature in the aging regime is the oc- 
curence of triangle relations between correlations at dif- 
ferent time scales. These triangle relations have an ultra- 
metric structure analagous to that found in RSB. At suf- 
ficiently large times ti and it is possible to express the 
correlation at intermediate times, i2 < ^ < ^i, in terms 
of a function /, which depends only .on the correlations 
and has no explicit time dependence. c2l Mathematically 



C{hM)^ f[C{h,t),C{tM)] , 
and there is an inverse function / 



(86) 



(87) 



Cugliandolo and KurchanEH^H give a complete account 
of these triangle relations and the properties of the func- 
tions / and /. The function / can have fixed points a, 
such that a — f{a, a). Each of these fixed points consti- 
tutes a correlation scale. We make use of this approach 
to calculate the FDT violation factor X, in the long time 
dynamics of the quantum rotor model. 

The saddle point equations are written down to 0{u^) 
in Section [II. To obtain the appropriate equations in 



the long time regime, the response must be split into a 
stationary part and an aging piece, as in Equations ( [7^ ) 
and ([76|). At long times, the contributions to the dynam- 
ical equations from the bath kernels should be negligible, 
since the system has been in contact with the bath suf- 
ficiently long for one time quantitites to have reached a 
limit. What is studied here is how the two time quanti- 
ties evolve after the interaction with the bath has ceased. 
We assume WEB and WLTM, and treat integrals in the 
manner explained by Cugliandolo and Lozano.Ej Hence 
the equations we wish to solve are (for ti ^ ^2) 



Ygdtl 



= - — i^^(ii,t2) + Aii?(ti,i2) 



+2J^ / dtR{ti,t)R{t,t2) 
Jq 

+X4 / dtC{ti,t)^R{ti,t)R{t,t2) 
Jo 

--Xih? / dtR{ti,tfR{t,t2), 
-J Jo 



and 



C{h,t2) + \lC{h , t2) + K2Ritl,t2) 



+ lx2C{h,t2f - ):\2C{tut2)R{hMf 



-HiC{ti,t2fR{tx,t2) - -K4fl^Riti,t2f 



/ dt[C{h,t)R{t2,t) + R{h,t)C{t,t2)] 
Jo 

+ -A4 / dtC{h,tfR{tM) 
-J Jo 



+\i / dtC{h,tYR{ti,t)C{t,t2) 
Jo 

1 1"^ 
--Xih^ I dt [3R{ti,tfR{t2,t)C{ti,t) 
3 Jo 

+Rih,t)^C{t,t2)] , 



(89) 



where we drop the aging subscript (since the only parts 
of the correlation and response we work with here are the 
aging pieces) . The values for the coefficients are displayed 
in Appendix |^. 

The equations derived here are very similar to those 
that have been written down for the soft spin SK 
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model.SEi If we assume that at sufficiently long times 
only terms that are invariant under time reparametriza- 
tion,s i,ri}m}iiji [yaore explanation is given in Section 
, then the dynamical equations take 



VD 



the form 

= XiR{h,t2) + X2C{h,t2fR{ti,t2) 



dtR{ti,t)R{t,t2) 



+X4 / dtC{ti,tyR{ti,t)R{t,t2) 
Jo 

0^ XiC{h,t2) + lx2C{tut2f 



(90) 



-2J' 



dt [C{h,t)R{t2, t) + R{ti,t)C{t, t2)] 



+ -A4 / dtC{h,tyRit2,t) 
•J Jo 

/•oo 

+A4 / dtC{ti,tfR{ti,t)C{t,t2). (91) 
Jq 

None of the terms with coefficients proportional to 
enter into these asymptotic equations, hence the depen- 
dence on quantum effects is only through the coefficients 
Ai, A2, and A4. There are also classical terms (in the 
sense that they do not have a coefficient proportional to 
h^) that are not reparametrization invariant. 

The equations above have extra terms relative] to the 
terms considered by Cugliandolo and Kurchan,c3 which 
are those with coeffcient A4 - we can recover their results 
when these terms are unimportant. The reason that we 
have these terms is that we include all terms of 0{u^) 
before looking at the long time limit. Following Cuglian- 
dolo and Kurchan we can convert the equations into a 
manifestly time reparameterization invariant form by in- 
troducing functionals F[C] and H[C] 



F[C] = - f\c'X[C'] , 
Jc 

H[C] = - [ dC'C'^X[C'], 
Jc 



(92) 
(93) 



and postulating a modified version of the FDT to relate 
the correlation and response 



Rih,t2) = i3x[cih,t2)]^c{h,t2)e{h - 12) 

OT2 

^p-^F[C{hMmh-t2) ■ 

Ot2 



(94) 



The relation ( |94| ) is also invariant under time reparam- 
eterization. We follow a similar reasoning to Ref. 40 to 
show that within an ultrametric scale 



X[al] 



X2al {2J^ + X4q^)2 
P (2J2 + A4af)i 



(95) 



We also find that within an ultrametic scale X[C] — 
X[al], and that the scales have vanishing measure and 



there is thus a continous ultrametric solution at correla- 
tion scale al , provided the following condition is satisfied 

2a:iF[a*] - X2 - X2af 



X2al + (3X4F[al]{2alX[al] - F[al]) 



> 0, 



where 



X2^Xi+ 2X4(3 I dC C'F[C'], 



I 2.P + A4g2 
2 J2 -I- X4af 



(96) 

(97) 
(98) 

(99) 



The results here ipeduce to those found earlier for the 
classical SK modeled in the limit that A4 ^ 0. Hence the 
solution for the FDT violation factor X[C] is 



X[C] 



X2C (2j2-t-A4g')2 
(i (2j2 + A4C2)i 



(100) 



Note that if we define q to be the value of the correla- 
tion for which X\C] = 1, then this implies q = /3/A2 = 



D. Reparametrization symmetry 

In deriving Equations ( pO| ) and ( ^ ) we assumed that 
the equations at long tinies are reparametrization invari- 
ant. As shown recently,E3 this assumption is justified, 
since at long times the dynamical equations governing 
the system flow to a "fixed point" of the reparametriza- 
tion group at which the equations are invariant un- 
der time reparametrization transformations. The time 
reparametrizations are transformations of the form 



t^i= h{t), 



(101) 



where h(t) is a differentiable function with dh/dt > 1 
(so as to stretch time). Transforming time also leads to 
transformation of the correlation functions - for a two- 
time correlation function G{ti,t2), we define a transfor- 
mation G ^ G such that 



dh 



^) G[h{t,),h{t2)], (102) 



where Af 2 are the scaling dimensions of the correlator G 
under the rescaling of the time co-ordinates ti^2- We term 
Af and A2 the advanced and retarded scaling dimen- 
sions respectively. The "fixed point" dynamical equa- 
tions are a set of equations which are invariant under an 
RpG transformation, i.e. if they are satisfied by G, then 
they are also satisfied by G. The idea of irrelevancy un- 
der KpG transformations, introduced in Ref. where 
certain terms in the long time limit scale to progressively 
negligible perturbations about an RpG fixed point, is the 
means used to obtain the dynamical equations at long 
times. 
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VI. DISCUSSION 

We have studied spin glasses of quantum rotors from 
a purely dynamical perspective. In the dynamical ap- 
proach, the phase transition between the paramagnetic 
and glassy phase can be identified by looking at the crit- 
ical slowing down of the dynamics. Alternatively, we 
also identify the transition at the breakpoint where it is 
no longer possible to satisfy the QFDT relation between 
correlators and response functions in the solution of the 
equations of motion. We find that, if one insists on a TTI 
solution, it is necessary to include a singular {S function) 
piece in the correlation to balance the fluctuation dissi- 
pation relation. The amplitude of this singular term is 
proportional to q, and this simple minded TTI solution 
is the analog of the usual replica symmetric equilibrium 
solution. 

For the glassy phase, we study the non-equilibrium dy- 
namics in the aging regime, where TTI is broken. We 
find that the relation between correlation and response 
for the quantum rotor model has the same character as 
in the classical case - specifically the OEFDR has the 
same form, but different coefficients from the classical 
case. We show that all of the terms in the dynamical 
equations that are only present in the quantum version 
of the model (with an explicit h dependence) are not in- 
variant under time reparametrizations, but instead they 
are RpG irrelevant and do not contribute in the long-time 
limit. The quantum terms, however, alter the short-time 
dynamics and consequently renormalize the coefficients 
of the classical (RpG invariant) terms in the long-time 
dynamical equations. 

These findings help us understand rather more pre- 
cisely how the behaviour of a spin glass system is modi^ 
fied by the introduction of quantum dynamics. BhattEiJ 
has argued that the scalings at the spin glass transi- 
tion should be the same as in classical models, due to 
the different timescales involved - f3h is the timescale 
for the quantum case (this corresponds to a frequency 
oi ujq — 8 X lO^-'^THz, where T is the temperature in 
Kelvin) , which at finite temperature will be much shorter 
than other relevant time scales in the problem. Another 
way to argue this is to look at the QFDT, which in the low 
frequency (long time) limit corresponds to the classical 
FDT. However, the glassy systems we consider here are 
never at equilibrium, and hence the long time limit in the 
non-equilibrium system is less obvious. One suggestion is 
that the QFDT should be generalized in a manner similar 
to the generalization to an OEFDR that differs from the 
classical FDT by introducing an effective tempecature 
(3X as the argument of the hyperbolic tangent .Ej i^he 
picture we discuss here and in recent previous workta is 
that at long times, due to reparametrization invariance 
of the dynamical equations, the dynamics become clas- 
sical in nature (as we would have naively expected from 
equilibrium); however, the effective temperature is differ- 
ent to that obtained in the corresponding classical model. 



due to quantum effects renormalizing coefficients. 

In essence, we find that in models where quantum dy- 
namics is introduced in a model with trivial spin commu- 
tation relations such as the rotors (in contrast to SU (2) 
spins, for example) quantum effects are only important 
in basically two aspects of the problem. The first is an 
equilibrium issue: the redefinition of the phase bound- 
ary due to the addition of quantum fluctuations, which 
further disfavors an ordered phase. The second has a 
non-equilibrium character: the renormalization of the 
classical terms in the aging dynamical equations, which 
translates into a renormalized effective temperature in 
the OEFDR. In this sense, this class of quantum spin 
glasses only correspond to "renormalized" classical sys- 
tems. The situation for truly SU (2) quantum spins, how- 
ever, is not clear and deserves further work. 

In addition to the quantum-classical connection that 
we partly address in this paper, there are also interesting 
issues connected to the time reparametrization invariance 
of the long time dynamics. The time reparametrization of 
the classical equations of motion appears to be a generic 
feature of classical spin glass models. The ipfpiianGe, p, 
the SK model was noticed by several authors ,c3'E3e3"E3'E3 
and we naturally encounter it in the rotor model as well. 
The form of the OEDFR seems to be a generic conse- 
quence pi the retarded and advanced scaling dimensions 
of C, i?J13 which in turn may be fixed by the form of the 
FDT. In this sense, it would be interesting to discover 
whether it is possible to construct different KpG fixed 
points, with different scaling dimensions. 
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APPENDIX A: TTI SOLUTIONS FOR THE 
RESPONSE 

One can attempt a solution for the response of the 
form R{uj) = P{u}) + iQ{uj), It is then possible to solve 
for P and Q and determine the quantum critical point, 
by writing x{ll}) = a{uj) + ib{Lu). The solution is 

R{uj) = + j5) ± ^v/(a^ ~ b'^ ~ 8J2) + 2iab, (Al) 

which can be separated into real and imaginary parts: 

X cos I tan-1 [ , | ) , (A2) 
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Qiiu) 



X sin tan 



2ab 



Va2 - 62 _ SJ2 



(A3) 



We obtain the critical point when the square roots in the 
arguments of the sine and cosine vanish. 



APPENDIX B: COEFFICIENTS FOR THE 
AGING EQUATIONS 

The coefficients that appear for the dynamical equa- 
tions in the aging regime are summarized below: 

Ai + u{l + M)q + 2(1 + M)u^B^ + 4J^Xoo 

+ ^(1 + M)uW^ - 1(1 + M)u-^'^Z (Bl) 
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